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^ . Abstract. In this short review paper, we discuss some of the recent advances in 

the field of parton fragmentation processes into hadrons as well as their possible 
modifications in QCD media. Hadron production data in e^e~, deep inelastic 
scattering and hadronic collisions are presented, together with global analyses of 
fragmentation functions into light and heavy hadrons and developments on parton 
' fH ' , fragmentation in perturbative QCD at small momentum fraction. Motivated by 

the recent RHIC data indicating a significant suppression of large-p^ hadron 
I ' production in heavy-ion collisions, several recent attempts to model medium- 

modified fragmentation, e.g. by solving "medium" evolution equations or through 
Monte Carlo studies, have been proposed and are discussed in detail. Finally we 
mention the possibility to extract medium-modified fragmentation functions using 
photon-hadron correlations. 



> 

, 1 Introduction and summary 

ON 

The fragmentation process - that is the transition from a highly virtual time-like parton, 
^ ^QCD' ^ coUimated bunch of hadrons - is a subject of first importance in QCD. 
Similarly to structure functions in deep-inelastic scattering, the momentum distribution of 
hadrons produced in the final state cannot be determined from perturbation theory. The reason 
is twofold: 

— the collinear splitting of one parton into two is divergent; 

— hadronization is a soft, non-perturbative, process. 



; I ■ Because of this soft and collinear sensitivity, the fragmentation function (FF) of a given parton i 

' with virtuality Q into hadrons h carrying an energy fraction x, D^{x, Q^), cannot be determined 

from first principles and should be extracted from experimental measurements. However, as long 
as Q is large as compared to j^qcdj the collinear divergences can be systematically factorized 
to all orders [1], 



Dlix, Q2) = ^ dz K^iz, Q\ Ql) X Dlixlz, Ql), (1) 

where the evolution of the parton i into the parton fc, with virtuality Qo and energy fraction 
z, is described by the evolution kernel K^{z,Q'^,Qq) and can be computed in perturbative 
QCD (pQCD). Therefore, even though the theoretical objects D^{x,Q'^) are strictly speaking 
infinite, their scaling violation is under perturbative control and has been successfully tested in 
hadron production in e^e'^ collisions, e.g. at LEP [2]. 

Based on these e+e^ data, several parametrizations of fragmentation functions into light 
hadrons (BKK [3], BFGW [4], KKP [5], Kr(etzer) [6], AKK05 [7]) and into photons (BFG [8]) 



^ Laboratoire d'Annecy-le-Vieux de Physique Theorique, UMR5108 



2 



Will be inserted by the editor 



have been performed. Despite the large amount of e+e" data, those FF sets still suffered from 
rather large uncertainties, especially at large x and/or small Q^. The situation has however 
considerably improved recently: experimentally, new measurements in deep inelastic scattering 
(DIS) and hadronic collisions at RHIC have been made available, while on the theoretical side 
two of the most recent fits, HKNS [9] and DSS [10], now include error analyses (work is in 
progress for AKK08 [11]), in analogy to what has been performed for parton distribution func- 
tions (see e.g. [12-14]). In the heavy-quark sector, the impressive results from the B-factories 
experiments (BaBar, Belle, CLEO) arc expected to bring in the near future significant con- 
straints on the non-perturbative component of FF into heavy hadrons. 

At small values of x ^ 1, or cquivalently large values of ^ = ln(l/a;), the Modified Leading 
Logarithmic Approximation (MLLA) [15-17] proved very successful in describing a variety of jet 
observables, such as single inclusive spectra or multiplicity distributions (for a review on MLLA 
phenomenology, see [18]). The success of this approximation, dealing with perturbative quark 
and gluon degrees of freedom, supports the idea of Local Parton Hadron Duality (LPHD) [19] 
which basically assumes that the shape of parton spectra remains virtually unmodified at the 
hadronization stage (only an overall normalization factor, K, = 0(1), is necessary to accom- 
modate the hadron data). In some cases though, typically whenever hadrons are produced at 
large angles with respect to the jet axis, the MLLA predictions fail to describe the experimen- 
tal measurements. Interestingly, the inclusion of terms of order O (a^) in the solution of the 
evolution equations, i.e. strictly beyond MLLA, appears to improve significantly the agreement 
with data [20]. 

Fragmentation processes have also received a lot of attention in the context of high-energy 
heavy-ion collisions where a dense and hot partonic system may be produced. Indeed, unlike 
in p-p collisions, the quarks and gluons produced with large transverse momenta in the initial 
nucleon-nucleon binary collisions in heavy-ion scattering (i.e. when the two incoming nuclei 
overlap) propagate through the strong colour field whose size and lifetime are 0(10 fm) 
produced from the soft "underlying" event. Therefore, the full fragmentation process of hard 
partons should be affected by the presence of the dense QCD medium, at least in principle. 
A typical example is parton energy loss caused by medium-induced multiple gluon emission, 
a process which has been widely studied over the past decade (see e.g. [21] for reviews). This 
mechanism may be responsible for the spectacular "jet quenching" phenomenon reported by 
the PHENIX [22,23] and STAR [24,25] experiments at RHIC, that is the suppression of single 
inclusive hadron spectra by a factor of 5 in central Au-Au collisions with respect to (properly 
scaled) p-p scattering. A lot of effort is now devoted in understanding how the parton multiple 
scattering process may affect QCD evolution. Based on the extensive phenomenology developed 
for jets produced in the "vacuum", i.e. in e+e~, DIS and hadronic collisions, several attempts 
to model fragmentation functions in heavy-ion collisions have been siiggested, cither by solving 
medium-modified evolution equations [26-30] or through Monte-Carlo studies [31-33]. 

In this short review paper, we will highlight the recent advances on fragmentation functions^ 

as well as their possible modifications in QCD media (see also [35] for a comprehensive overview 
of that subject). Section 2 is devoted to fragmentation processes in the vacuum. After discussing 
the new data available (Section 2.1), emphasis is then put on the recent FF sets and their 
improvements with respect to older studies (Section 2.2), the developments beyond MLLA 
at small x (Section 2.3) as well as parton fragmentation into heavy hadrons (Section 2.4). 
Moving to medium-modified fragmentation processes in Section 3, data in electron nucleus and 
nucleus-nucleus collisions are first summarized (Section 3.1). After having briefly recalled the 
basics of parton energy loss processes (Section 3.2), we discuss the modelling of QCD evolution 
in presence of a QCD medium (Sections 3.3 and 3.4) and review recent attempts to develop 
parton showers in nucleus- nucleus collision in Section 3.5. Finally, we show in Section 3.6 
how performing photon-hadron momentmn correlations may help to extract medium-modified 
fragmentation functions in heavy-ion collisions. 



^ The field of spin-dependent fragmentation functions and transversity is not addressed in the present 
paper. For a review on that subject, see Ref. [34]. 



Will be inserted by the editor 



3 



2 Fragmentation function studies 
2.1 Data 

2.1.1 e+e^ annihilation 

The measurements of the total fragmentation function in e+e^ annihilation, or scaled-energy 
distributions, 

F'\x, s) = — ^(e+e- ^ hX) (2) 
cTtot dx 

of light hadrons were first performed at DORIS/PETRA at DESY (ARGUS [36], CELLO [37], 
JADE [38], TASSO [39] experiments) at centre-of-mass energies ^/s — 12-36 GeV, as well as 
at SLAG PEP (HRS [40], MARK II [41], TPG [42]) at = 29 GeV, and at KEK TRISTAN 
(TOPAZ [43]) at y/s = 58 GeV. More recently, hadron production in e~^e~ collisions has also 
been measured at the Z° pole {y/s = 91.2 GeV) both at GERN LEP-I (ALEPH [44-48], 
DELPHI [49,50], OPAL [51-53]) and at the SLAG linear coUider (SLD [54,55]), and at higher 
energy, = 133-209 GeV, at LEP-II (DELPHI [56,57], L3 [58], OPAL [59,60]). 
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Fig. 1. BaBar preliminary measurements of scaled-energy distributions of data vr*. A"*, and p/p in 
e^e~ collisions at — 10.54 GeV. Taken from Ref. [61]. 

The distributions have been measured for both unidentified and identified (tt^, K^, K^, 
p/p, A/ A) hadrons on a wide range of momentum fractions, from a; ~ 5 10^"^ {£, ~ ^ lux ~ 5) 
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at LEP/SLC up to a; ~ 0.8, the statistical uncertainty becoming however fairly large at high 
X. Hadron spectra inside gluon jets have boon obtained from a selection of symmetric thrcc-jet 
events, e+e~ qqg and using fe-tagging techniques [45,47,50,62,63]. Distributions in giuon 
jets arc found to be softer than the ones from quarks jets (also reported are stronger scaling 
violations [63]) because of the larger colour charge of the gluon: the higher multiplicity making 
hadrons to be produced at smaller x by momentum conservation. Moreover, several experiments 
(DELPHI [49], OPAL [53], SLD [55]) were able to perform flavour-tagged measurements and to 
disentangle light-quark (uds) from c-quark and 6-quark jets; note in particular that the OPAL 
experiment reported on hadron spectra in each individual light flavours [53] . A discrepancy was 
found between the NLO calculations using the FF parametrizations then available (BFGW, 
KKP, Kr) and the bottom and gluon jet data, despite a good agreement with flavour-inclusive 
and udsc-ict measurements [63]. For more detail about these data samples, we refer the reader 
to the global fit analyses [3-11] and the references therein. 

Impressive preliminary measurements on hadron spectra in e+e" annihilation have also 

been reported very recently from B-factories (BaBar, BELLE, CLEO) around the T(4S)-mass, 
^/s = 10.58 GeV, and below {y/s < 7 GcV) through initial-state radiation studies, e~^e~ 
7 /i X. While those experiments are naturally suited to study heavy-quark hadrons (discussed 
in Section 2.4), it is worth to mention that light-hadron spectra can also be measured with an 
accuracy better than what has been achieved so far. For the sake of an illustration, preliminary 
measurements of tt*, , and p/p inclusive spectra by the BaBar collaboration [61] are plotted 
in Figure 1 together with older data from ARGUS, TASSO, and SLD. 



2.1.2 DIS 



Inclusive hadron spectra have also been studied experimentally in DIS events, mostly in e p 
collisions at HERA by HI [64,65], HERMES [66], ZEUS [67,68] experiments but also in neutrino 
scattering, i/jj^p and V^p, by NOMAD [69]. The virtuality range probed at HERA, 10 < Q < 
100 GcV, coincides fairly well with e+e~ data, while being somewhat smaller (1 ^ Q ^ 10 GeV) 
with NOMAD. 

In order to ease the comparison with e+e~ collisions, the DIS events are analyzed in the 

Breit frame for which the photon has a virtuality but no energy. In the parton model, the 
struck (anti)quark has longitudinal momentum pz = Q/2 and Pz ~ ^Q/2 before and after 
the collision with the virtual photon. The; hemisphere containing particles with negative Pz is 
called the current hemisphere and should be directly comparable to e^e^ collisions in which a 
quark is produced with energy \/s/2, while the target hemisphere {pz > 0) describes the proton 
remnants. 

Quite remarkably, the momentum distribution (2) of inclusive charged hadrons has been 
measured recently by HI for 9 different intervals in x and for 6 bins in virtuality Q, allowing 
for a systematic study of scaling violations of fragmentation functions [65]. In all the (x, Q^)- 
bins for which these data can be compared to e+e~ results, a very good agreement is obtained 
supporting the universality of the fragmentation processes (see also Refs. [2,70] for universality 
checks from e+e~ to DIS) expected from factorization theorems. A small difference between 
the two scattering systems occurs at low energy scale, coming most probably from next-to- 
Icading order corrections present in DIS and absent in e+e~ annihilation. Comparing data with 
theoretical expectations reveals that NLO QCD predictions, using either KKP, Kr, or AKK05 
FF sets, are not able to describe properly the Q-dependence of HI data, the scaling violation 
being much stronger in the experiment [65] . 

Another successful check of universality was given by the recent high-precision ZEUS mea- 
surements of multiplicity distributions, dA^'/drich at different scales [68] , which compare very well 
with TASSO and LEP data. In particular, these data confirm that the product {rich) x diV/drich 
is a function of rich/ (rich) independently of the jet energy, a property expected in QCD and 
known as Koba-Nielsen-Olesen (KNO) scaling [71], as can be seen in Figure 2. 
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Fig. 2. KNO scaling of charged hadron multiplicity distributions measured in DIS events by ZEUS 
and compared to e^e~ data taken by TASSO and LEP experiments. Taken from Ref. [68]. 



2.1.3 Hadroproduction 

Because of the structure of both the projectile and the target, inclusive momentum hadron 
distributions in hadronic collisions may a "priori reveal little on fragmentation processes, since 
for instance no scaled-momentum distributions such as Eq. (2) can be constructed. However, it 
happens that the recent high- precision measurements of the BRAHMS [72], PHENIX [73,74], 
and STAR [75-77] collaborations in p-p collisions at RHIC, and CDF [78,79] in p-p colhsions 
at the Tcvatron, actually allow for additional constraints to be set on FF parametrizations (see 
discussion in the next Section 2.1.4). 

The transverse momentum spectra of pions, kaons and (anti)protons at mid-rapidity^ were 
measured in those experiments with a high accuracy up to roughly ~ 10 GcV, and even up to 
■p^ ~ 20 GeV for neutral pions by the PHENIX collaboration [74] . Perhaps even more interesting 
are the measurements of strange particles such as K^^ [72,75] and Aj A [77]. Finally, the RHIC 
experiments BRAHMS and STAR were able to determine spectra differences between positive 
and negative-charged pions, zAtt^ = d(T(pp tt+X) — Aaipip —> 7r~X), a quantity obviously 
vanishing in e+e" annihilation and in p p collisions at mid-rapidity. Because of the wealth 
of new measurements available from RHIC, those data sets are now systematically included 
in the recent global fit analyses DSS [10] and AKK08 [11] (but not HKNS [9]) discussed in 
Section 2.2. Before these data, the inclusive particle spectra in hadron collisions taken e.g. at 
the CERN SppS by the UAl [80] and UA2 [81] collaborations were only used to test a posteriori 
the fragmentation functions extracted from e+e~ data [5,6], with the exception of the BFGW 
set [4] which was actually the first analysis to consider hadroproduction measurements in order 
to get better constraints on the gluon fragmentation functions. 

Single inclusive charged hadron spectra inside jets have also been measured very precisely 
in p-p collisions at the Tevatron by the CDF experiment [82] on a large variety of dijet masses, 
ranging from 80 to 600 GeV. Finally, let us note that jets of energies up to 50 GeV have been 
reconstructed lately by the STAR collaboration at RHIC [83]. This will eventually allow for 

^ Measurements at large positive rapidity, y ^ 3, are also reported by BRAHMS [72] and STAR [76]. 
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single momentum spectra to be measured down to rather small values of x and for different 
hadron species which was not done at the Tevatron. Preliminary measurements have recently 
been shown [84] which shapes seem to agree qualitatively well with the MLLA expectations 
discussed below. 

2.1.4 Comparing systems 

After having discussed the various experimental measurements related to fragmentation pro- 
cesses in e^e~ , e{u)p, pp, and pp collisions, it is worth to discuss the respective advantages and 
drawbacks of each collision system. 

Electron-positron annihilation is the cleanest process and seems by far the most suited 
reaction to study fragmentation functions. Indeed, to leading order in the strong coupling 
as, the momentum distribution F''(x,s), Eq. (2), measured in e^e^ is simply related to the 
individual parton-to-hadron fragmentation functions, D^ix, Q^), entering the fixed-order QCD 
calculation: [85] 

Ff^{x,s)= J2 9i{s) D^{x,s). 

i=q.q 

where gi{s) is the electroweak coupling. In analogy with the space- like case, the observable 
F^{x, s) is reminiscent of the structure function ^2(0;, Q^) measured in DIS, and the individual 

FF D'j' (x,.s) the analog of parton distribution functions, f^{x,Q'^). Beyond the leading order, 
factorization theorems [1] ensure that F^^ can be expressed as 

F\x,s)=Y. f ^cA^,s,iA DUz^I^'), (3) 

Jx Z Z 

up to higher- twist terms which are power suppressed, e.g. O ^^qcd/'^)- Both the parton-to- 
hadron FF, D'^, and the so-called coefficient functions, Ci, are factorization-scheme dependent 
(e.g. MS or DIS) and depend on an arbitrary fragmentation scale, /i, taken to be O (v^) in order 
to avoid large logarithmic corrections, ln(/i^/s). The coefficient functions have been computed 
in the MS scheme a long time ago at NLO [86] and more recently at NNLO [87] . 

On top of its simplicity another obvious asset of e+e" scattering is the wealth of data 
available, as discussed in Section 2.1.1. Moreover, the fragmentation of heavy quarks into either 
light-quark or heavy-quark hadrons is not as suppressed (at least as long as mg <C Vi) as e.g. 
in DIS processes. These various reasons explain easily why almost only e+e~ data have been 
used, until recently, to constrain FF parametrizations. 

Despite those attractive features, however, e+e^ collisions constrain pretty weakly the giuon 
FF which appear at NLO, as can be seen in Eq. (3). Also, since hadrons produced in e+e~ come 
from the fragmentation of both the quark and the anti-quark jet (at LO), those data only allow 

for the extraction of flavour-inclusive FF, e.g. + (= Dq^'^ from charge conjugation). 
This is clearly at variance with DIS and p-p collisions in which quarks are produced more 
abundantly than anti-quarks in the hard process because of the valence-quark distribution in 
the (projectile and) target hadron - at least as long as the momentum-fraction at which PDF 
are probed is not too small, say xbj > 0.1. Therefore, the individual and shall be 
determined independently from DIS and p-p measurements^. 

As already discussed in Section 2.1.2, the advantage of DIS over hadroproduction exper- 
iments is the close; (;onnection offered with e+e^ scattering, which allows one to observe the 
universality of fragmentation processes. Also, fragmentation functions can be probed on a wide 
(a;, Q^)-range unlike in e+e~ for which the hard scale is set at LO by the centre-of-mass energy 
of the collision, Q = y/s/2. Note however that at NLO, the hard scale also depends on the 
typical angle between two neighbouring jets 9, Q = O {Ejg^sm{9/2)) [15,88] which allows for 
the variation of the hardness Q in 3-jet events in e+e~ collisions [89]. 



^ In principle, rapidity asymmetries in p-p collisions could also be used. 
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Coming to liadroproduction experiments, p-p and p-p scattering are ideal in order to con- 
strain the gluon fragmentation functions, especially when xbj -C 1 hence when the hadron p^^ 
is small as compared to Note also that FF extracted in hadronic collisions are typically 
probed at large values of x (typically x ^ 0.7 0.8 at RHIC [90]) because the parton densities fall 
steeply when xbj = O {pj_/x) is large. Therefore the kinematic region probed in hadronic col- 
lisions is different, and hence complementary, to the one in e+e~ collisions. Finally, measuring 
hadrons inside jets makes it possible to compare scaled-momentum spectra in hadroproduction 
experiments together with e^e^ and DIS and check further the universality of fragmentation 
(see e.g. Figure 4 below). 



2.2 Global fit analyses 



The many attempts to extract FF from e+e" data have already been mentioned in the Intro- 
duction. We would like in this Section to focus only on the three most recent parametrizations, 
namely HKNS [9], DSS [10], and AKK08 [11] which therefore should be used when computing 
inclusive hadron production in elementary collisions emphasizing in particular the different 
assumptions made in each global fit analysis. For clarity, the assumptions and ingredients of 
each set are summarized in Table 1. 



Table 1. Various characteristics of recent global fits analyses AKK08 [11], DSS [10] (and DSV [10] 
for strange particles), and HKNS [9]. For each set we list the data samples used, the strange particles 
available, the error analyses, and further theoretical details. 



AKK08 



DSS/DSV 



HKNS 



e'^e / pp / pp e+e / e p / pp 
K°,A + A K'^,A + A 

— Lagrange errors 
mass effects — 

large-a; resum. — 



Hessian errors 



2.2.1 Data samples 



As far as data samples are concerned, all three sets used almost the same e+e" measurements 
taken at LEP/SLC as well as at lower energy. Therefore, differences arise rather on the inclusion 
or not of DIS and hadroproduction data. While HKNS considered e+e~ only, both AKK08 and 
DSS included the RHIC results (BRAHMS, PHENIX, STAR) to get a better handle on gluon 
FF; AKK08 also taking into account measurements of and A/ A spectra performed at ^/s = 
630 GeV by CDF [79] . DSS is the only analysis which included semi-inclusive DIS measurements 
by HERMES [66], despite the rather low Q"^ ~ 1-5 GeV^ probed in that experiment. 



2.2.2 QCD evolution 



The DGLAP equations [91-93] governing the Q"^ evolution of parton-to-hadron fragmentation 
functions read 

dDi{x,Q'^) _ 



dlnQ2 



^Pn{-^,c,sm) Dj{z,Q') (4) 
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where Pji are the NLO time-like splitting functions determined in [94,95]. These evolution 
equations are either solved through a "brute- force" z-integration (HKNS) of (4), or using the 
standard Mellin techniques (AKK08, DSS) writing Eq. (4) as 

j 

in momentum space where the moment 

' dz 



D^{N,Q')= I ^ z''Diiz,Q'). 



Jo 

can then be Mellin- inverted to obtain the fragmentation function Di{x,Q'^), 

where C is a contour which lie at the right of all singularities. The reader may refer to [96] for 
a discussion on the various approaches. Mellin techniques also proved useful in computing the 
hadroproduction cross sections at NLO in a very efficient way. Indeed, the Mellin-transform of 
the hadronic coefficient functions C^ j{z), i.e. complicated functions involving the product of 
two parton densities and the partonic QCD cross sections, need only to be computed once [10]. 
Therefore the hadron-hadron production cross section, symbolically written as 



i,j,k ^ u C i A 



can easily be determined, exploiting the fact that the moments of the FF drop rapidly with Af, 
D%{N) oc A''-^ [10], where /3 is the large-x power law of the FF, oc (1 - zY . 



2.2.3 Fit assumptions 



All groups assume SU(2) isospin symmetry for the sea (also called unfavoured) fragmentation 
functions D'l, that is processes for which the hadron h does not contain the flavour i as a 
valence parton. Taking the 7r+ = \ud) as an example, they assume that 

HKNS assume moreover that Df = Dl \ yet the SU(3) isospin symmetry could be pretty 
badly broken because of the strange quark mass; DSS report for instance = O.SSiP^ from 
their fit. In the valence (favoured) sector, HKNS and AKK08 conjecture that SU(2) symmetry 
also holds, — DJ , while DSS allow for the different normalization of the valence plus sea 
quark FF, D^_^^ = A/ x Du+u-, A/" being a free parameter which is not necessarily equal to one 
(as assumed by AKK08 and HKNS). Prom their global fit analysis, and in particular thanks 
to the SIDIS HERMES preliminary measurements, DSS obtained M = 1.10, yet the value of 
A/" depends somehow on the individual fitted data samples. They also check the assumption 
made by Kretzer [6] that unfavoured fragmentation functions are suppressed at large x by 
one extra-power in (1 — x) with respect to favoured ones, e.g. Dg^(x) = (1 — x)D^^ (x). This 
working hypothesis turns out to be in rather good agreement with 7r+ data from HERMES [66] , 
although it somehow failed to describe tt" production, (x) 7^ (1 — x)D^ {x) [10]. 
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The analytic parametrization of the FF at an initial scale Qo {1 GeV^ for DSS and HKNS, 
2 GeV^ for AKK08) is 



Dt{x,Ql) = Ntx^^\l-xf' 



(5) 



in AKK08 and DSS, and taking the simpler form = in HKNS. The normalization param- 
eters iV/* are constrained in such a way to satisfy the following momentum sum rule, 



E 

h 



Jo . 



(6) 



for all parton species i and all Q^. Rather than giving the FF into charge-sign identified hadrons, 
e.g. ^ and D'^ , AKK08 chooses another "basis" , the inclusive-charged DJ '^'^ as well as the 
"valence" quark FF, = D'^^ — , which is needed to compute charge-sign asymmetries 
e.g. at RHIC or in semi-inclusive DIS. 




Fig. 3. HKNS parametrization for charged pions, Df /2, at NLO for all flavours, compared to the 
AKK05, KKP, and Kr sets previously available. The bands indicate the theoretical uncertainty of the 
HKNS fragmentation functions. Taken from Ref. [9]. 
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2.2.4 Error analyses 

Besides the use of new data sets available from RHIC for instance, the most important im- 
provement of the recent FF sets concerns the treatment of errors in the global fit analyses. 
Fragmentation studies have in particular greatly benefited from the progress made in the ex- 
traction of parton distribution functions over the last few years [12-14]. 

Using the Hessian method [12] which assumes a parabolic dependence of the function for 
parameter values away from its minimum, HKNS compute the spread of the FF propagating the 
(correlated) errors on the individual parameters entering Eq. (5). The advantage of this error 
analysis is to make rapidly clear which fragmentation functions are the most uncertain and 
in which kinematical domain^. For the illustration, the HKNS individual NLO fragmentation 
functions into charged pion data are displayed in Figure 3 and compared to the "older" sets 
AKK05, KKP, and Kr; the envelopes of the FF indicate the theoretical error from the Hessian 
analysis. 

DSS used instead the Lagrange multiplier technique [13] which helps to understand the 
range of variation of physical observables allowed from the global fit, taking the truncated 
second moment as an example in [10]. Also, interestingly, the influence of a specific data sample 
on a given observable can be easily determined within this approach. 

2.2.5 Hadron mass corrections and large-a; resummation 

Unlike DSS and HKNS, AKK08 have not performed an error analysis, although work towards 
this goal is in progress. What is performed in AKK08 and not carried out elsewhere, however, 
is the treatment of hadron mass eflFects, rrih ^ 0, in the extraction of fragmentation functions. 
They are taken into account by making explicit the kinematical differences between the light- 
cone scaling variable x entering the fragmentation functions and the energy (xe) or momentum 
(xp) fraction used experimentally e.g. in e+e" collisions (corrections can be taken care of in a 
similar manner in hadronic collisions): 



The AKK08 fits have been performed for all hadron species letting ruh as a free parameter. 
Remarkably, the fitted mass turned out to be amazingly close to the true masses (1% difference 
is quoted) for baryons {p/p, A/ A), yet 10% larger for pions, an observation attributed in [11] 
to the decay of the p-meson, whose mass is much larger than that of the pion; note however 
that the fitted kaon mass proves 30% smaller than its actual value. On the phenomenological 
side, hadron mass effects are visible at low x and/or small s/s. 

The other novelty of the AKK08 analysis is the large- iV (large x) resummation of leading 
logarithms, In""*"^ N, and next-to-leading logarithms, a" In" A^, of the Mellin transform of the 
NLO e+e~ coefficient functions, Ci, appearing in Eq. (3) [97]. Large-a; logarithms appearing in 
the DGLAP evolution are also resummed to LL and NLL accuracy, using [98] . Even though the 
inclusion of large-x resummation improves the fit for , p/p, and especially A/ A production 
(x^/ndf = 1.45 instead of x^/ndf = 1.73 without resummation), differences with the purely 
fixed-order analysis are however pretty small, and only visible either at low energy or above 
x > 0.7. Given that fragmentation functions are extracted within a resummed analysis, one 
may also wonder whether AKK08 should in principle be used in a fixed-order calculation. 

2.2.6 Comparing FF sets 

A detailed comparison of the various sets would not be appropriate here; the reader is rather 
referred to the original publications of each analysis [9-11] for a complete discussion. The 

* Note also that a FORTRAN code released by HKNS allows for computing the theoretical uncertainty 
coming from their FF set on any observable. 
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general features are thus only briefly sketched in this Section. First of all, given the amount 
of e^e~ data available, the fragmentation of light quarks into tt^ is rather similar in each set. 
However, the gluon fragmentation functions, Dg , reported by AKK05/AKK08 and especially 
that of DSS (both using hadronic data) are much harder than the one given by HKNS which is 
only extracted from e+e~ measurements. This can be seen in the upper left panel of Figure 3 
where AKK05 and HKNS are compared. Therefore, potentially largo differences in the QCD 
predictions of hadron -spectra e.g. at the LHC could be expected, AKK08 and DSS being 
a priori more trustful than HKNS because of the additional constraints brought by RHIC 
and Tevatron data. Nevertheless, the slight disagreement between AKK08 and DSS found at 
high a: is a hint that systematic theoretical uncertainties on the g tt* procc^ss may still be 
large. Indeed, even though the hadronic measurements help to constrain gluon fragmentation 
processes, the data in p-p and p p reactions arc not as precise as the e^e~ data samples, and 
uncertainties remain substantial at high values of x. Even though the agreement between data 
and theory may vary from one set or one hadron species to another, all NLO predictions of 
TT^ and spectra reproduce fairly well the available data. Remarkably, those sets prove also 
rather successful in predicting baryon production such as p/p (despite the agreement is not 
as good as for tt^), yet the FF into p/p are rather dissimilar in each set. In the case of A/ A 
production, however, a large discrepancy between theory and STAR data in p-p collisions [77] 
is apparent either when using AKK08 [11] or DSV [10] sets, which both undershoot the data 
by roughly a factor of 5. Finally, the largest differences between the parametrizations are seen 
for the charge-asymmetry performed at RHIC which should eventually be able, when the data 
become more precise, to discriminate the various assumptions regarding favoured vs. unfavoured 
fragmentation functions. 

While it is a common procedure to vary the renormalization/factorization in order to give 
an estimate of unknown higher order corrections, that is the perturbative error, it is a clear 
that the predictions of hadron spectra should ideally also use all the three recent sets AKK08, 
DSS, and HKNS, in order to quantify the "non-perturbative uncertainty" of these calculations. 
Having a tool similar to the recent LHAPDF interface for parton densities [99] would of course 
help significantly in this respect. 



2.2.7 Towards NNLO analyses? 



All the global fit analyses performed recently showed that the agreement between the FF sets 
and the current data was (much) higher when going from LO to NLO accuracy. Therefore it is 
legitimate to wonder whether even higher-order corrections would improve further the situation 
or not. 

The calculation of the space-like splitting functions at three loops, O (a^), has been carried 
out a few years ago by Moch, Vogt and Vermaseren [100], then allowing for the first analyses of 
parton densities at NNLO accuracy [101]. Unfortunately, no complete calculation is yet available 
at this order for time-like splitting functions. At leading order, MS space-like and time-like 
evolution kernels are identical, an observation known as Gribov-Lipatov reciprocity [91,102]. 
Beyond that order time- like splitting functions can be determined from an analytic continuation 
in X [95, 103] or through the conjecture made by Dokshitzer, Marchesini and Salam [104] who 
suggest a simple correspondence between space-like and time-like evolution to all orders from 
a redefinition of the standard DGLAP equations''. 

Based on those two proposals, Mitov, Moch and Vogt recently achieved the computation of 
the NNLO non-singlet splitting functions [105]. There is therefore a real hope that analyses for 
fragmentation functions at NNLO can be performed in the near future. 



^ Basically, the scale at which the PDF or FF must be evaluated in the r.h.s. of (4) is Q^z and /z, 
respectively, with new "universal" splitting functions V{z,as{Q^ /z)). 
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2.3 Fragmentation at small x 

The DGLAP evolution (4) deals with large collinear logarithms which are resummed within the 

Leading Logarithmic Approximation (LLA) . As a consequence of strong transverse momentum 
{k±_) ordering in the cvohition, all powers of hard collinear logarithms, as In fc^ with k±/u) 
1, are resummed. This scheme describes the space- like and time-like evolution of the parton 
densities and fragmentation functions, respectively, which both satisiy DGLAP equations. In 
the literature it is also called the "fixed order approach". However, when radiated gluons in 
the process reach the infrared region (that is, .x <C 1), the rcsummation of large soft logarithms 
becomes also necessary, which leads to the Double Logarithmic Approximation (DLA) [15,106]. 

In DLA, as a consequence of strong angular ordering (AO) in cascading processes, all powers 
of soft and collinear logarithms, ln(l/.T) Infc^ arc resummed. This scheme describes the 
time-like evolution of fragmentation functions and takes properly into account the crucial QCD 
coherence effects: soft and collinear gluons interfere destructively once k± becomes smaller than 
a certain value. In addition, it was demonstrated that the dominant contribution to inclusive 
particle production is obtained when the successive emissions of soft and collinear gluons are 
strongly ordered in their energies (uji ^ Wi+i) and emission angles {"di 3> This hierarchy 

in energies and angles in DLA gives rise to a simple probabilistic interpretation in terms of 
classical parton shower cascades, which can easily be implemented in Monte Carlo approaches. 
Indeed, within this picture, branching processes in a high energy jet can be displayed in terms 
of Feynman diagrams at tree level, which tremendously simplifies the problem by getting rid of 
all possible interferences. However, in this approximation the energy balance at each vertex is 
ignored (basically assuming 1 — a; ~ 1), which leads to an overestimate of particle production 
in DLA. 

Adding the double logarithms (DLA) to the single collinear (but hard) leading logarithms 
(LLA) leads to the Modified Leading Logarithmic Approximation (MLLA) proposed in Refs. [15- 
17]. As compared to DLA, MLLA partially restores the energy balance for each gluon emission - 
although it remains approximate - and takes explicitly into account the running of the coupling 
constant in the evolution, from the hardest scale Q to threshold Qo- Remarkably, the proba- 
bilistic interpretation encountered in DLA is preserved in MLLA with a simple prescription on 
the angular ordering, which now becomes exact: 'di > "ffi+i, after azimuthally averaging over 
the phase space of the emitted gluon. Consequently, all powers of ag ln(l/a;) Ini? + ag Int? are 
resummed in the MLLA, where ln(l/.T) Ini? = 0(1) and a^ln?? = O (y^oJ). The collinear 
terms are known to partially restore the energy balance that is overestimated in the DLA. 

The small-x evolution equation can be obtained from the DGLAP evolution equation with 
an appropriate change of the evolution variable, from the virtual mass-squared to Q^, 

dD,{x,Q^) 



din 02 



(i.i'Q'). P) 



which reduces to Eq. (4) at x = O (1). Because of this change of scale in (7), the anomalous 
dimensions jji are no longer divergent at j = 1 (which corresponds to a; = 0). In the gg case, 
it is given by [85] , 

^-^^■'"^^=^^i-l + 4,(,,a.)' 
with cig = asNc/2-K. Performing the Taylor expansion of (8) in (j — l)^/as 4; 1, [85] 

lgg{j,as) = v^- i(j - 1) + ^j={j -If + O {^^^jPj , (9) 

allows one to extract the leading small-a; contribution to multiplicities or particle spectra. Note 

that; interestingly, 7gg oc ^/oi^ at small- .t. On the contrary "/gg can be written as a power series 
in as when Taylor expanding (8) at small coupling constant (and/or large j — 1), 

7s9(i,a.) = aa j^-a«'(j^Y)3+^((7^)" 
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10 100 

Mjj sin(ec) GeV/c2 

Fig. 4. Location of the peak of the single-inclusive distribution xD{x, Q^) as a function of the energy 
scale Q, measured in e^e~ , DIS, and p-p collisions. The prediction in the Leading-Logarithmic Ap- 
proximation (short-dashed line), MLLA (solid line), and neglecting coherence effects (long-dashed line) 
are shown for comparison. Figure taken from Ref. [82]. 



From the anomalous dimension (9) the moments of the fragmentation function, D(j,t) ~ 
can be obtained [85], 



2tt 



77 5/2 



(11) 



where is evaluated at the scale and b — (llA^c/3 — 2n//3)/27r. Performing the inverse 
Mellin transform of Eq. (11) leads to the DLA fragmentation function. 



xD{x, Q^) ~ exp 



2a- 



(12) 



which is Gaussian in ^ = ln(l/2:). The distribution - known as "hump-backed plateau" - has 
a maximum because of gluon QCD coherence at small x (large ^), see Figure 11 below. The 
location of the peak and the width of the distribution are given by 



(13) 
(14) 



The DLA inclusive spectrum (12) is correct at asymptotic energies. When including single 
logarithm corrections, the MLLA spectrum is softened as compared to DLA and its peak is 
shifted towards larger values of ^ by an amount oc y/a^, [15] 




= In 



Q 



A, 



1 

2+« 



32iVc7r 



with a = llNc/'i + 2nf/3N^. The MLLA spectrum is expressed in [15] in terms of the inverse 
Mellin transform of a confluent hypergeometric function. Recently, an integral representation 
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has been obtained by Perez Ramos beyond the limiting spectrum approximation, that is when 
the infrared scale Qq ^ ^qcd = 1ii(Qo/^qcd) 0)i [10'''] 

tt+,+A)//i^..-. r^x (w^) M-"*' 

J J (27ri) Jo J^ + s \[LU + s)iyJ K'^ + -V 

(with y = \n{Q / A^^-^) — ^) which was then estimated analytically in Ref. [108] using the steepest 
descent method. 




12 12 3 



ln(k^/1GeV) ln(k^/1GeV) 

Fig. 5. CDF preliminary measurements [110] of fcx-distributions for hadrons inside a jet compared to 
MLLA (dashed line) and NMLLA (solid line) calculations. Figure taken from Ref. [20]. 

The shape of the MLLA inclusive spectrum describes successfully the momentum distribu- 
tions of hadrons measured in e+e^, DIS and p-p collisions at the Tevatron (see e.g. [18] for a 
comprehensive review on MLLA phenomenology). In particular, the decrease of parton spectra 
at small x (large ^) because of QCD coherence is observed experimentally. Another check of 
the MLLA framework is the scale dependence of the location of the hump - i.e. independent 
of the spectrum normalization - which compares well with the experimental results (Figure 4) 
and brings further evidence for coherence effects in QCD at small x. 

Remarkably, apart from the value of ^qc^ j the only parameter entering the MLLA spectrum 
is an overall normalization factor K. of the order of 1. The similarity of the shape of partonic 
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spectra computed in MLLA on the one hand and that of hadrons measured experimentally on 
the other hand supports the idea of Local Parton Hadron Duality (LPHD) introduced in [17,19]: 
hadronization happens locally and thus does not distort dramatically the spectra of quarks and 
gluons when they turn into hadrons. 

More differential observables have been considered lately, such as the transverse momen- 
tum spectra of hadrons inside a jet, dN/dk±, which has been computed for the first time 
by Perez Ramos and Machet at MLLA and in the limiting spectrum approximation {Qq = 
^qcd) [109]. More recently, the calculation has been extended beyond the limiting spectrum 
and including O (ag) corrections which are strictly speaking next-to-MLLA (NMLLA) [20]. 
On the experimental side, the CDF collaboration reported on preliminary measurements of 
/c_L-distributions of inclusive-charged hadrons from dijet production in p-p collisions at a/s = 
1.96 TeV [110], on a large variety of scales, Q = M^j sin(6'c/2) = 19 155 GcV, where M^j is the 
dijet invariant mass and 9c the jet cone radius. As can be seen in Figure 5, the NMLLA predic- 
tions reproduce very well the experimental measurements (both data and theory are normalized 
at ln(A:_L/l GeV) = —0.1) for all energy scales. Another interesting fact is that MLLA spectra 
prove harder than CDF data, which is an indication that recoil effects neglected in MLLA are 
necessary to account for the experimental results. 

Recently, the small-x evolution equations in the coherent branching formalism [111] (of 
which Eq. (7) is an approximation) has been solved numerically [112], allowing for the calcula- 
tion of single spectra with exact energy conservation in the evolution. Although no significant 
differences with MLLA are observed at large jet energies (ln((5/(5o) ^ 7), the shape of the 
hump-backed plateau proves somehow modified at lower energy [112]. 

Finally, another interesting development has been suggested by Albino, Kniehl, Kramer 
and Ochs [113] in order to unify the standard large- a; DGLAP evolution equation, Eq. (4), 
together with the DLA evolution at small values of a;, Eq. (7). The DGLAP equation Eq. (4) 
is solved using new splitting functions which resum the double logarithms to all orders (DLA 
part) supplemented by the regular part P of the fixed-order splitting function in order to avoid 
any double counting: [113] 

as P\z)^P^'^{z,as)+asT^{z). 

Remarkably, solving DGLAP equation with these modified splitting functions allows for an 
excellent description of single inclusive e^e~ data at all values of x, from very small x (^ ~ 5) 
and up to very high values, x ~ 0.9. Moreover, the extrapolation from large to small x does 
not require the help of LPHD and the inclusion of an arbitrary normalization factor. 



2.4 Fragmentation into heavy-quark hadrons 

Heavy quark production can be computed perturbatively because the large mass, mg ^ ^qod) 
acts as a coUinear cut-off for the gluon emission by a heavy quark. Therefore, the fragmentation 

function, Of {x,Qo,mQ), of a parton with flavour i into a heavy quark and at an initial scale 
Qo can be computed order by order in perturbation theory. The calculation has been performed 
&tO{as) in [114]: 



D^{x,QlmQ)=5{l-x) + 



27r 



l-x 



log 



Ql 



(1 - xfm^Q 



D^{x, Ql mg) = ^ [x^ + (1 - xf] log ^ 



D'^^^{x,QlmQ) = 0{a1) 



(15) 

(16) 
(17) 



with Cp = 4/3, Tp = 1/2. The Dq fragmentation function receives a singular contribution 
corresponding to no gluon emission, while only appears at order as from the g — > QQ 
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splitting process f md = O (a^) from the double cascade q{q) q{q)g followed hy g ^ QQ. 
The initial scale is Qo = C {^q) in order to avoid large logarithms appearing in Eqs. (15) 
and (16). Once these initial conditions are known, the FF can be evolved through LO or NLO 
DGLAP evolution, Eq. (4), from Qq to the hard scale Q of the reaction. This procedure ensures 
that the leading a"(Q) In" Q/Qo and next-to-leading a^" (/u) ln"~^ Q/Qo logarithms of coUinear 
origin are resumed to all orders. 

As can be seen from Eq. (15), the initial condition for the heavy-quark fragmentation fimc- 
tion contains large terms, log(l — x)/{l — x) which are divergent in the x ^ 1 limit. These 
terms arise at any order in the perturbative expansion, because of the non-cancellation of real 
and virtual terms in gluon emission, and need to be resumed in order to obtain an accurate 
description of fragmentation in this kinematic limit. This task has been achieved first at leading 
logarithmic [114] and then at next-to-leading order accuracy in [97,115], through the inclusion 
of Sudakov form-factors in the initial condition (15). 

As in the case of light hadron production, the heavy-quark scaled-momentum spectrum is 
given by 

FQ{x,s,mQ) = Y^ f ^CA-,S,A Df{z,t,^mQ). (18) 

^ Jx ^ ' 

Even though the sum runs over all flavours z, it is worth to mention that in e+e~ collisions 
heavy quark production from gluon splitting only contributes at a few percent level as indicated 
for charm production by the ALEPH data [116] (see also [117] for theoretical studies), even 
though it is expected to be somewhat larger in hadronic collisions where gluons are produced at 
the leading order. Because of the large mass, the heavy quark carries most of the jet momentum 
and therefore the single inclusive spectra of heavy hadrons are peaked towards large values of 
X [118], unlike the spectra of light hadrons that we have seen in Section 2.1 are peaked at small 

X. 

When computing heavy-flavoured hadrons, H, however, non-perturbative effects come into 
play because of the Q ^ H hadronization mechanism. Therefore the heavy-quark spectrum 
(18) has to be convoluted with a non-perturbative fragmentation function, i?np, m order to 
give the heavy-hadron spectrum: 

F"{x,s)=Y. [ T^'^G'''"''^) ^^^^ 
which in Mellin space simply reads 

F"{N,s) = Y, Ci{N,s,ij')Df{N,i^\mQ)D^p{N). 

i 

Several analytic parametrizations of non-perturbative fragmentation functions have been given 
(summarized in Table 2) which parameters can be constrained from a comparison with existing 

data. 

Similarly to the measurements of light hadron production discussed in Section 2.1, the 
spectra of bottom mesons have been measured in e+e~ collisions at the Z° pole by the LEP 
experiments [127-131] and by the SLD collaboration [132, 133] with an excellent statistical 
accuracy. In particular, both the OPAL and SLD data tend to exclude specific parametrizations 
of fragmentation functions such as Peterson et al., CoUins-Spiller and BCFY, favouring instead 
the Lund, Bowler, and to a lesser extent Kartvelishvili models [129,133]. Bottom production 
has also been measured in DIS by the HI collaboration [134] as well as in hadronic collisions 
by UAl [135] and more recently by CDF [136] and D0 [137] at the Tevatron. 

The measurement of charmed mesons in e+e" collisions at LEP energy is slightly more 
delicate than bottom production since momentum spectra need to be corrected for important 
B-meson weak decays. Still, rather precise measurements have been performed at LEP [138]. 
However the bulk of data regarding _D-meson production has been taken at lower energy. For 
quite some time, the only data available at y/s = 10 GeV by the ARGUS [139] and CLEO [140] 
collaborations had large statistical uncertainty and rather poor momentum resolution. Since 
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Table 2. Noii-pcrturbative heavy-quark fragmentation functions, L'np(a;). The fi appearing in the 
BCFY parametrization are polynomials of degree i. 



Model Parametrization 

BCFY [119,120] [i!'d::r.i^ [3 + Eti(-^074r)] 

Bowler [121] a;-(i+'''''' '"i>(l - a;)°exp {-bml/x) 

Colangelo-Nason [122] x^{l - x)" 

Collins-Spiller [123] 

Kartvehshvili [124] a:"''(l - x) 

Lund [125] a;"^ (1 - a;)" exp (—6 Tnj_/x) 

Peterson et al. [126] x~^ (l - a:"^ - (1 - x)~^)~^ 



then, high-precision measurements of the D"*", D*"'", and Z?*" spectra have been performed 
by BaBar [141], BELLE [142] and CLEO [143]. Finally, charm measurements were also carried 
out in DIS by HI [144] and ZEUS [145], yet with a slightly lesser accuracy. Interestingly, the 
D*"^ momentum spectra from ZEUS also disfavour the Peterson non-perturbative fragmentation 
functions [145]. 

The importance to get stringent constraints from e"^e~ data on the non-perturbative frag- 
mentation function (and in particular its first moments) to predict charmed and bottomed 
meson production in hadronic collisions has been highlighted recently in p p collisions at the 
Tcvatron. The CDF collaboration reported that their measurement of spectra was a factor 
of 3 in excess over the "QCD background" [146]. However, as shown in [147], this conclu- 
sion was actually drawn on the basis of a QCD calculation using the Peterson et al. FF that 
is excluded from LEP data. Predictions using more accurate non-perturbative fragmentation 
functions which reproduce well the first < 10 moments of the S-mcson spectra at LEP, 
somewhat increase the bottom cross section thereby reducing the excess of CDF data® to a 
factor 1.7. Such an enhancement could then easily be explained by the present uncertainties on 
parton densities and the scale dependence of the FONLL calculation [147]. 

In order to illustrate the phenomenology of heavy-flavoured production, we briefly sum- 
marize in what follows the main results of a recent and detailed analysis of D and S-meson 
production in e+e" collisions and the B-factories and LEP energy [148]. The analysis is per- 
formed using the NLO initial conditions (15)-(17) for heavy-quark FF together with NLO 
coefficient functions Ci, both terms being properly resummed to NLL accuracy to account 
for soft gluon emission. The NLO evolution of fragmentation functions is carried out with 
a proper matching conditions when crossing the bottom threshold [149], but neglecting any 
gluon splittings, g QQ, which are pretty much suppressed as already mentioned. Finally, the 
non-perturbative fragmentation function is taken from an average of a singular contribution, 
(X S{1 — x) and the Colangelo Nason parametrization [122] (see Table 2). Within this approach 
(and taking into account for possible initial-state radiation in e+e" collisions), the spectra of 
D*+, and D*+ D"/D*" prove in excellent agreement with BELLE [142] and CLEO [143] 
data at ^/s ~ 10.6 GeV. Interestingly, however, the extrapolation of these calculations to LEP 
energies turns out to be in contradiction with ALEPH measurements (which are the most pre- 
cise so far) at large x. The disagreement between theory and the ALEPH data is illustrated in 
Figure 6 where the data/theory ratio is plotted as a function of the moment of the fragmenta- 
tion functions . The reason for this disagreement may come from large power corrections to 



The final CDF measurements also slightly decreased as compared to the first preliminary results. 
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Fig. 6. ALEPH D*^ data compared to QCD predictions in momentum space, fitted by a simple 
functional form 1/(1 + 0.044(iV - 1)). Taken from Ref. [148]. 



the e+e~ coefficient functions, Ci/^/s or C2/S, whose strength can be fitted to ALEPH data, 
Ci = 0.57 ± 0.16 GeV or a somewhat larger value C2 = 5.1 ± 0.3 GeV^. Unfortunately the 
lack of data at intermediate scales between the T(4S') and the mass does not allow for a 
clarification of the scale-dependence of such corrections. 



3 Medium modifications 



After having reviewed the experimental and theoretical studies of fragmentation functions, we 
address in the second part of the paper the possible modifications of fragmentation processes 
induced by the presence of dense QCD media, such as nuclear matter in electron-nucleus scat- 
tering or quark-gluon plasma (QGP) in heavy-ion collisions. 



3.1 Data 



As we have seen in Section 2, the wealth of data in e+e^, DIS and hadronic collisions allow 
for detailed studies of fragmentation processes "in the vacuum". Unfortunately, the clean e^e~ 
process is obviously of no use in order to investigate possible medium modifications to frag- 
mentation functions. Therefore, one has to rely either on DIS experiments on nuclear targets, 
e-A, or heavy-ion collisions, A-A; the former being ideal to probe cold nuclear matter'' while 
the latter is devoted to the study of hot and dense - though expanding - QCD medium. Data 
samples involving nuclei in the initial state are however much more scarce and kinematically 
limited than those involving protons. 



3.1.1 DIS on nuclear targets 



The HERMES experiment at DESY has performed extensive measurements of hadron produc- 
tion in semi-inclusive DIS (SIDIS) on various nuclear targets: N, Ne, Cu, Kr and Xe [150-152]. 
Further, HERMES was able to perform the identification of 7r+, 7r~, tt", , K~, p, and p 
separately. The hadron yield in a given (z/, x)-bin, where v is the virtual photon energy in the 



Proton-nucleus scattering is also interesting to probe fragmentation in nuclear matter, yet it is not 
as clean as DIS on nuclei. 
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target rest frame and x = Eh/v the hadron energy fraction, is normalized to the number of 
DIS inclusive events, Nj^. Assuming collinear factorization to be valid for hadron production in 
SIDIS on nuclei, the leading-order ratio can be written as: 

^ i=g, q 1=9) <1 

(20) 

where Q is the virtuality of the photon and a"^ i = ''{x^^,!/) the LO 7*9 cross section. 

The fragmentation function d'^^^ indicates that the hadron is produced off a nuclear target, 
eA ^ e h X. When x is not too small, which is the case of the HERMES experiment, x > 0.02, 
sea-quarks in the target are small and therefore only up valence quarks contribute to the cross 
section (the contribution from down quarks being suppressed because of their smaller electric 
charge squared, = e^/4). In this approximation, the ratio simply reduces to the up-quark 
fragmentation function: 

,,,, 

It is indeed particularly interesting since the large current uncertainties on the nuclear parton 
densities, f^, appearing in (20) cancel out in this ratio^. The double production ratio in a 
heavy nucleus over that in a deuterium target ^, 



Nl diydx I 7V|, dvAx 



(22) 



(with Q"^ = 2mArXg. , tun being the nucleon mass) is then determined as a function of u or x. 
As an example, the quenching of charged hadron, pion and kaon spectra measured as a function 
of V by HERMES [152] are shown in Figure 7. 

Those data therefore prove ideal in order to study the detailed kinematic dependence of the 
medium- modified fragmentation process. This is will be of course only true if we can attribute 
the deviation of R\{z, v) from unity to be due to a modification of the up quark fragmentation 
into the hadrons. In that sense, the factorization formula Eq. (20) has to be seen as a model 
assumption. Suppose for instance that the typical time to form a hadron is of the order or less 
than the nuclear size, then inelastic rescattering of that very hadron in the nuclear medium 
may be responsible for the relative suppression of hadron production in e A scattering, i.e. 
R\{z, v) <1. Various models based on this "hadron absorption" mechanism have been proposed 
to account for the trend of the data [156-158]. 

The CLAS experiment at JLab is currently performing similar measurements [159], yet 
presently at somewhat smaller energies, v = 2-3 GeV, ^ 1 GeV^, for which perturbation 
theory may no longer be under control. Hopefully the upgrade at 12 GeV will allow for a 
more systematic study at higher energies {v = 2-9 GeV, = 2-9 GeV^) of all hadron 
species [159], including all light mesons and strange baryons, shedding new light on medium- 
modified fragmentation in cold nuclear matter. 



3.1.2 Heavy-ion collisions 

Large-pj^ hadron spectra have been measured in Au-Au collisions at RHIC (-^s^n = 200 GeV) 
by the BRAHMS, PHENIX, PHOBOS and STAR experiments. What is often reported exper- 
imentally is the ratio of the hadron yield at a given in Au-Au collisions (either minimum 

^ The first measurements of hadron attenuation in DIS on nuclei is due to the EMC collabora- 
tion [153]. Since hadron production was not normalized to the number of DIS events, the reported 
suppression in e-A with respect to e-p is actually the consequence of the modifications of parton 
densities, the so-called . . . "EMC efi^ect" . 

" We write D^-^ = in (22) since we do not expect any modification of fragmentation in such a 
small nucleus. 
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bias or in a centrality class C) over that in p-p collisions, properly scaled by the number of 
nucleon-nucleon (NN) binary collisions in that centrality class: 

The number ( A'^coii) | ^ is determined from the Glauber multiple scattering theory (see e.g. Ap- 
pendix I of [1601). It is expected to be (iVcon)l = 779 for the 20% most central Au-Au 

lc<20% 

collisions at RHIC energy. Since hard processes, such as large-p^^ production, are expected to 
scale with (A^coii)) the quenching factor (23) is R^{p^) = 1 in the absence of any "nuclear" ef- 
fect, which may come either from the properties of the initial state, e.g. modifications of parton 
densities in the incoming nuclei, either from final state interaction in a dense QCD medium. 

In central Au~Au collisions at y^Vn = 200 GeV, a significant suppression of neutral pions, 
R""" ~ 0.2 up to ~ 20 GeV has been reported by PHENIX and PHOBOS [162], a similar 
attenuation for ?7-mcsons up to p^ ~ 10 GeV being also reported by PHENIX [163]. STAR [25] 
(and also BRAHMS [164] and PHENIX [165]) observes a similar quenching of the inclusive 
charged hadron yields {p^ < 10 GeV). High-p^ baryons {p, p, A, A) appear also suppressed, 
but at p^ values higher than for mesons [166]. Interestingly, the hadron suppression pattern 
is less pronounced when going to more peripheral collisions [25]; in particular, i?'' < 1 for the 
most peripheral (C = 80-92%) collisions, see Figure 8. Similar features are also observed at 
y/s = 130 GcV [22,24,167]. For a complete discussion of RHIC measurements in heavy-ion 
collisions, we refer you to the comprehensive reviews by BRAHMS [168], PHENIX [169] and 
STAR [170], or to Ref. [171] for a more concise survey of large p^^ measurements at RHIC. 

These experimental results, sometimes referred to as "jet quenching" [sic], certainly one 
of the most spectacular phenomenon observed at RHIC, is often quoted as one big piece of 
evidence for parton energy loss in QGP discussed in the next Sections. What is particularly 
remarkable indeed is the absence of suppression in d-Au collisions where the quenching factor is 
around unity [172]. This is the indication that the mechanism responsible for the pion quenching 
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Fig. 8. Quenching factor of neutral pion production in Au-Au collisions at = 200 GeV for 

various centralities measured by PHENIX. Taken from Ref. [161]. 



is due to final state interactions - hence for which parton energy loss is a natural candidate - 
and cannot be attributed to e.g. strong modifications of the parton densities in the colliding 
nuclei. 

What happens then at lower energy? The limited phase-space available (^s^n — 20 GeV) 
at the SPS fixed-target facility does not allow for a systematic study of "large" hadron 
spectra. The WA98 7r° measurements indicated an enhancement of the yield in Pb-Pb central 
collisions [173]. However, a re-analysis of the 7r° reference spectra in p-p scattering leads to a 
quenching factor ~ 1 or even below 1 in the most central collisions [174]. Hopefully the 
data at RHIC intermediate energies will clarify this issue [175]. 

To study parton energy loss and medium-modified fragmentation functions would ideally 
require the reconstruction of jets in heavy-ion collisions; the huge background makes of course 
this task highly delicate. Nevertheless, thanks in particular to important theoretical develop- 
ments on the jet reconstruction algorithms in a high-multiplicity environment [176], preliminary 
measurements by STAR [177], but also studies at the LHC by ALICE [178], ATLAS [179], and 
CMS [180], look very promising. 



3.2 Parton energy loss 

3.2.1 Induced gluon spectrum 

Many signatures for quark-gluon plasma (QGP) formation have been suggested. Among them, 
it was argued first by Bjorken in the early eighties [181] that hard quarks and gluons propagating 
through QGP may experience multiple scattering and therefore lose some energy because of 
the induced gluon emission. Experimentally, the immediate consequence would be a depletion 
of high-energy jets in heavy-ion collisions with respect to p-p scattering. 

This process was then revived a decade later when Thoma and Gyulassy [182] and Gyulassy, 
Pliimer, and Wang [183] computed perturbatively the radiative energy loss of high-energy 
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partons in QGP. Since then various approaches have been developed to determine the gluon 
radiation spectrum, dJ/dw, of hard partons undergoing multiple scattering. 

Baier, Dokshitzer, Mueller, Peigne, and SchifiF (BDMPS) developed a perturbative frame- 
work to describe the medium-induced gluon emission process from the soft multiple scattering of 
hard partons in cold [184] and hot [185] QCD matter. The calculation assumes that the number 
of collisions, or opacity, is large: n = L/X » 1, where L is the medium length and A the parton 
mean free path. At small energy the typical time to produce gluons exceeds its mean free 
path, t ^ uj^^ > A. Consequently several scattering centres act coherently to stimulate gluon 
emission, giving rise to the Landau-Pomeranchuk-Migdal (LPM) effect in QCD. Because of the 
destructive interference, the LPM spectrum is suppressed in the infrared, wdZ/dw oc 
as compared to the independent Bethe-Heitler gluon spectrum, uidl/du) oc ui~^. The BDMPS 
approach is equivalent [186] to the powerful path-integral formulation of Zakharov [187], later 
generalized by Wiedemann [188] for an arbitrary number of scatterings. Another approach by 
Gyulassy, Levai and Vitev has been carried out for thin media, n = 1,2..., from which the 
gluon spectrum at any opacity is obtained recursively [189]. Let us also mention the finite- 
temperature approach by Arnold, Moore and Yaffe (AMY) [190, 191] for moderately "hard" 
partons with energy O (T), which was first developed to address the thermal photon rate in hot 
QGP [191]. Finally, a twist expansion has been proposed in Ref. [192]. Since the goal in this 
paper is to address mostly fragmentation processes and their modifications, we refer the reader 
to Ref. [21] for a more comprehensive discussion on the diff'erent energy loss formalisms. 



3.2.2 Quenching weights 



In order to model the effect of energy loss on parton fragmentation into hadrons, an essential tool 
is the probability distribution in the energy loss, or quenching weight. Neglecting interference 
eff'ects among the radiated gluons, which are 0{as), Baier et al. [193] gives the quenching 
weight a simple Poisson expression: 



oo n 
n=0 ■ Li=l 



dlicoi) 



i=l 



exp 



/ 



dto 4— 
do; 



(24) 



in which the gluon spectrum is the only ingredient entering the calculation. The quenching 
weight helps to make the connection between the parton energy loss in QCD on the one hand 
and the experimental consequences on the other hand. In the general case, 7^(e) has a discrete 
contribution on top of a continuous part p{e) [194] , 



V{e)=pa <5(e)+p(e). 

The probability for no-gluon emission is given by pq 
number of gluons radiated with an energy larger than oj, 



(25) 

lini exp[— iV(a;)], where N{ijj) is the 



I 



For asymptotically large media the BDMPS gluon multiplicity is infrared divergent, N[uj) oc 
w"^/^, and therefore the probability for no energy loss vanishes. With proper kinematic bounds 

for the gluon emission in finite-length media, however, the spectrum becomes finite in the 
infrared and therefore the quenching weight acquires an explicit singular contribution, po 7^ 
in (25) [194]. 

Using an integral representation of Eq. (24), the quenching weights have been computed us- 
ing the BDMPS and Wiedemann induced gluon spectrum [154, 194] and were made available for 

practical use either in terms of a mmicrical program [194] or an analytic parametrization [154]. 
For the sake of an illustration, the energy loss distribution is plotted in Figure 9 for various 
quark energies (in units of the energy loss scale, Wc)- 
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Fig. 9. Quenching weight using the BDMPS gluon spectrum for several quark energies in units of the 
energy-loss scale Uc- Taken from Ref. [154]. 

Recently, Peshier [195] proposed a new way to determine the elastic, or collisional, quenching 
weight based on (24). As a first order Markov process, the momentum distribution of partons 
produced at a time t -\- 5t is given by that at time i. 



in discrete momentum space, where the transition matrix Tqp is simply related to the p q 
transition rate. Therefore the probability for a hard parton with energy E{p) to lose an energy 
e = E{p)—E{q) in the time interval 5t is basically given by Tqp which, by successive iteration (i.e. 
matrix multiplications), gives the full quenching weight. On top of being numerically fast, one 
of the advantages of that method is that it should be appropriate for any parton energy, while 
the Poisson approximation (24) explicit breaks down when considering finite parton energies as 
discussed in [154]. 

3.2.3 Expanding media 

The medium-induced gluon spectrum, and therefore the quenching weight, is determined for 
static and uniform media. However, the dense QCD medium produced in the early times of a 
heavy-ion collision experiences fast (and mostly longitudinal) expansion. In order to account 
for this dynamical expansion, the transport coefficient characterizing the scattering power of 
the medium, q, is usually rescaled according to [196], 



where a characterizes the time-dependence of the medium energy-density, nir) cx t~", and 
q{T) = <?(ro) (to/t)". The purely longitudinal (or Bjorken) expansion corresponds to a = 1, and 
is often assumed in phenomenological applications. Recently, a remarkably simple prescription 



fq{t + St)=%p fp{t) 




(26) 
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has been given recently by Arnold in order to determine d//dcj in a finite and expanding 
medium [197]. Therefore, applying this recipe using e.g. hydrodynamical space-time evolution 
will allow eventually for the computation of more realistic quenching weights without relying 
on the assumption (26). 



3.3 Kinematic rescaling 

The first approach to model medium-modified fragmentation functions accounting for parton 
energy loss in QCD media is due to Wang, Huang and Sarcevic [198]. The model basically 
assumes that the evolution of the fragmentation process is not affected by the medium. 
The effect of energy loss is simply to rescale the initial parton energy, from say E to E — e, 
after which the parton fragments as if in the vacuum (see sketch of the model in Figure 10). 
Therefore, the medium-modified FF, call it T>, can be written as [198]: 

xV^{x,Q')^ deV{e)x* Dl{x*,Q') + deV{e) ^ Dl ^^,Q'y (27) 

where the energy- fraction carried away by the hadron is shifted from x to x* = x/{l — e/E^). 
Because FF are steeply falling functions of a; (see Section 2), even a small shift in the momentum 
fraction could lead to a significant hadron yield suppression. The second term in the r.h.s. of 
Eq. (27) accounts for the fragmentation of the radiated gluons induced by the medium. Although 
it satisfies the momentum sum rule (6), it is legitimate to question the use of the large- a; FF and 
at very low in order to fragment those gluons into hadrons (a simpler LPHD prescription 
may appear more appropriate). Further, a hard parton which has lost an energy e through the 

medium is likely to have emitted several gluons {n ^ 1) whose energies sum up to e, uji — e, 

i—l.n 

rather than a single gluon with energy uji — e (see e.g. [154]). 




Fig. 10. Schematic view of the energy loss process and the medium-modified fragmentation function 
in the energy rescaling model. 

This model - or variants of it - has been widely used to describe the hadron quenching 
measured at RHIC [199,200]. The trend of the data as a function of was in particular 
nicely reproduced when geometrical effects - partons being produced anywhere in the dense 
medium - are taken into account [200]. The calculation based on the BDMPS/ Wiedemann 
framework, however, led to values of the transport coefficient somewhat larger [200] than what 
is expected in perturbation theory [184,201] . More recently, progress has also been made towards 
the implementation of a full hydrodynamical expansion of the produced medium - constrained 
by soft/global observables - in energy-loss scenarios, first in [202] and more recently in [203] 
aiming at a consistent description of soft and hard probes in heavy-ion collisions. 

In order to illustrate the qualitative features of the model, let us take the case of SIDIS 
reaction for simplicity. Forgetting about the convolution for the sake of the discussion, assuming 
V{€) — 6{€ — (e)) as long as x is not too large, the suppression of hadrons using (27) in (22) 
becomes, 

R\ix,.)^l + '-^'-^^. 

1/ ax 
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First of all, the suppression is an increasing function of the photon, hence of the; up cjuark energy 
in the nuclear rest frame. In particular, note that in the Bjorken limit u ^ (e) , hence Q ^ (e) at 
finite x^^., the effects of the final state interactions become negligible, V^{x,Q^) ~ D'j'{x,Q^). 
Secondly, it is is interesting to notice that the suppression of hadrons is sensitive to the (log- 
arithmic) slope of the ordinary fragmentation functions [155]. Hence, since the unfavoured 
fragmentation functions are softer than the favoured ones (see discussion in Section 2.2.3), this 
in turn leads to a larger suppression of hadrons which do not contain up valence quarks (as long 
as x^. is not too small). This effect is actually seen in the HERMES data (see Figure 7) which 

report a stronger suppression of K~ than K~^, < , in agreement with the energy loss 
model calculations-'^^ [155]. Finally, at large x the phase-space to emit gluons with energy less 
than (1 — x)Ey becomes very much restricted. Therefore the suppression shall be all the more 
stronger as ,t gets larger, as can be seen from the upper limit of the integral (27). In the a: ^ 1 
limit, in particular, the DIS quenching factor is given by the probability for no-gluon radiation 
in the nuclear medium, 

R'Xix l,u)=po. 

The x dependence from the rescaling model turns out to be as well in rather good agreement 
with the HERMES measurements [155]; see also Ref. [158] for a more recent discussion. 



3.4 Parton evolution in QCD media 

The above-discussed rescaling model appears reasonable to describe the energy loss of the 
leading parton. However, the rescattering of the full parton shower in the medium is not properly 
taken into account in the Poisson approximation, Eq. (24), of the quenching weights. Moreover, 
the dependence of the parton energy loss process is not taken into account in that model (the 
medium-modified T>{x,Q^) (27) actually satisfies the standard DGLAP evolution equations). 

In order to address these issues, several approaches using the MLLA resummation scheme, 
encountered in Section 2.3, have been proposed recently [26-28]. The common idea of these 
studies is to modify the evolution kernel, Kj in Eq. (1), in order to model the extra gluon 
radiation in the case of radiative energy loss [26,27], or the rescattering of the emitted partons 
in the medium [28]. 

In the approach of Borghini and Wiedemann [26] , the gluon radiation induced by the medium 
modifies the singular part of the QCD splitting functions; for example. 



2(1 + /mod) 

— {i- + Z) 

(1-Z + 



(28) 



where /med is a parameter controlling the "amount" of induced gluon radiation. In particular, 
taking /med brings back to the ordinary splitting functions. The more infrared singular 
behaviour of the induced gluon spectrum, with respect to that in the vacimm, justifies the 
choice of modifying only the 1 /z part of the splitting functions. Shifting splitting functions in 
order to reproduce parton energy loss in QCD media was first advocated in [192]. 

The medium-modified MLLA equation can therefore be solved in a very same way as Eq. (7) 
in the vacuum in order to compute a variety of observables in presence of a dense QCD medium. 
The single inclusive spectrum D'^{x,Q^) has been determined in particular in [26] and is re- 
produced in Figure 11. As can be seen, the usual hump-backed plateau is "distorted" by the 
medium: the number of highly-energetic particles (small ^) is quenched as compared to the 
expected vaciuim spectrum, while the soft gluon yield is dramatically enhanced at large ^. This 
is nothing but the manifestation of the (approximate) momentum-conservation in the MLLA 
evolution: the medium converts part of the energy of a hard parton into the production of many 
soft partons (there is no "energy loss" as such!). This model was recently extended beyond the 
limiting spectrum approximation {Qq ^ ^qcd)j that is for hadrons with mass m ~ Qq [205]. 



^° The absorption models also predict a stronger suppression of K than because of a larger 
inelastic rescattering cross section in nuclear matter. See e.g. the short discussion in [204]. 
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Fig. 11. Single inclusive hadron spectrum as a function of ^ using vacuum splitting functions (compared 
to e^e~ data) and using medium-modified splitting functions. Taken from Ref. [26]. 

A specific dependence on the hadron species inside the jet ("jet hadrochemistry" ) is predicted 
and could be tested with ALICE at LHC which has good identification capabilities [178], as 
well as CMS [180] to a certain extent. 

In the approach of Armesto, Cunqueiro, Salgado and Xiang (ACSX) [27], the medium- 
modified splitting functions are directly related to the medium- induced gluon spectrum, d/™°'^/dz dQ^, 
computed in the Wiedemann parton energy loss framework [188], 



As compared to the Borghini- Wiedemann model, the medium modifications are now explicitly 
dependent of the parton virtuality through the induced gluon spectrum. No energy loss effects 
are expected at the time of hadronization which occurs on long time scales. Therefore, the FF 
at the initial scale Qo ^ 2-4 GeV^ is given by that in the vacuum (specifically, the KKP set [5] 
is assumed) [27]. A virtue of this approach is that at large virtualities, Q ^ ^qcd; rescaling 
model (27) using the quenching weights (24) is formally recovered [27]. 

The medium-modified FF have been determined on a large range of a; = 5 x 10~'^-1 and 
for all [27]. As expected, the effects of parton energy loss is to soften the fragmentation 
functions, with a large suppression at high- a; and a clear enhancement at small x, see Figure 12. 
Of course, the larger the medium length or the transport coefficient, the stronger the medium 
effects. What is perhaps less intuitive is the dependence: the effects of parton energy loss 
are more pronounced as gets larger-'^^. This behaviour is explained in [27] by the longer 
time-evolution, even though the time spent by the hard parton in the medium may naively be 
set by the medium size and lifetime. Within the same framework, the medium dependence of 
the mean and the dispersion of multiplicity distributions of partons inside quark and gluon jets 
has also been obtained analytically recently [29,30]. 

The recent RHIC data on the suppression of single inclusive electron spectra in Au-Au 
collisions [206] has triggered a renewed interest on elastic, or collisional, energy loss in QCD 

Note that medium-modified FF do not coincide with those in the vacuum for asymptotically large 
and at finite x as one could have naively thought for higher twist effects. 



P{z) 



as dzdQ^ 



(29) 
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Fig. 12. Medium-modified gluon-to-pion fragmentation functions for different values of q and medium 
lengths. Taken from Ref. [27]. 

media [207]. In the same spirit as the above describe radiative energy loss models, it was 
suggested in [28] to take into account 2^2 scattering processes by the inclusion of a scattering 
function, /C, in the leading order evolution kernel. The function /C is proportional to the medium 
gluon density and to the gg — > gg scattering cross section. 

At " ^^^ ' 4 H + m|,)2' 

where the Debye mass, mo, acts as an infrared regulator. A large suppression is predicted 
at high X but, unlike the radiative energy loss approach [27], no enhancement is expected at 
small values of x. For consistency, however, it would be interesting to include such a term, 
K, = O (a^), in a DGLAP equation at NLO, and to check whether the medium effects remain 
similar. 

Finally, let us also mention that fragmentation functions including parton multiple scattering 
have been determined in a higher-twist framework [192]. Although the calculation was first 
applied in the context of electron-nucleus collisions with phcnomcnological success, it has also 
been extended for the case of hot QCD media. Recently, Aurcnche, Zakharov, and Zaraket 
(AZZ) [208] however showed that the correct gluon emission spectrum actually vanishes when 
the number of rescattering is n = 1; this claim is disputed by Wang in [209], yet maintained by 
AZZ in [210]. 

3.5 Medium-modified parton showers 

Parton showers (PS) based on Monte Carlo techniques, e.g. HERWIG [211] or PYTHIA [212], 
share many advantages. First of all, they allow for the energy-momentum conservation through- 
out the full evolution, unlike analytic approximations such as DLA or MLLA. Moreover, PS 
are able to characterize the underlying event activity in the collisions - an important issue with 
heavy-ions - and to ease the comparison with experimental measurements. Finally, they also 
allow for a better understanding of the microscopic dynamics. 
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Fig. 13. n-jet distribution predicted by JEWEL with elastic (left) and radiative (right) energy loss. 
Taken from Ref. [31] 

In heavy-ion collisions, several parton shower models which describe the jet evolution in 
QCD media have been developed recently [31,32]. In [31], the Monte Carlo programme JEWEL 
(Jet Evolution With Energy Loss) is first tested for various observables, sucdi as jet shapes 
in e+e" collisions. Concerning the medium modifications, radiative energy loss is implemented 
using the model by Borghini and Wiedemann [26] described in Section 3.4 where /mod is a free 
parameter. Elastic rescattcring is also included in JEWEL through the aforementioned 2^2 
scattering process [28]. Therefore, exploring various jet observables with JEWEL may allow 
for both energy loss mechanisms to be disentangled. A good example is the n-jet fraction as 
a function of the resolution scale log 2^0(7/0) assumed in the jet algorithm. Figure 13 shows for 
instance that elastic rescattering almost has no effects on the vacuum n-jet distribution, while 
on the contrary the number of jets produced at a given yo increases somehow when radiative 
energy loss are taken into account [31]. This is an encouraging step towards the understanding 
of the mechanisms responsible for jet quenching. Another PS [32] was also developed from a 
modification of the PYTHIA shower algorithm PYSHOW [213]. The effect of parton energy 
loss is to increase the virtuality of the rescattering partons (thus making them to radiate more) 
by an amount AQ"^ which directly depends on the space-time dependence of the transport 
coefficient q = q(jx., t) given by a 3D hydrodynamical evolution. Finally, note that the approach 
of ACSX [27] is being implemented in PYTHIA as well and preliminary results have already 
been given [33]. 



3.6 Measuring (medium-modified) fragmentation functions 

There is an obvious need to identify observables in heavy-ion collisions which may be sensitive 
to parton energy loss in hot QCD media. As stressed e.g. in [200], the single inclusive hadron 
-spectra is unfortunately not to sensitive to the microscopic dynamic underlying the "jet 
quenching" phenomenon. As mentioned in Section 3.1, there is a real hope that measuring 
spectra of hadrons inside reconstructed jets, yet difficult, may be possible. In this Section 
another observable is explored: the double inclusive production of a prompt photon and a 
large-pj^ hadron. 
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Consider the leading-order QCD Compton process for instance, qg 57. Because of mo- 
mentum conservation, the photon is produced back-to-back to the hard quark, with equal and 
opposite transverse momentum. As a consequence, the momentum imbalance variable. 




reduces to the fragmentation variable, z^^ — x, in this leading-order (LO) kinematics. There- 
fore, there should be in principle a clear connection between the experimentally accessible 
photon-hadron momentum-imbalance distributions on the one hand, and the theoretical quark 
fragmentation function into hadrons on the other hand, as first suggested in [198]. 

There a a few caveats though. First of all, the photon can itself be produced by the coUinear 
fragmentation of a leading parton. In this case, there is no correlation between the hadron 
and the photon momenta since the two jets fragment independently. This is precisely the most 
important "background" channel which can get rid of using appropriate isolation criteria. Other 
limitations which could somehow complicate the picture are higher-order corrections as well as 
soft and coUinear gluon radiation. 




In order to illustrate the interest of this observable, the 7-/1 distributions computed in QCD 
at LO are shown in Figure 14 for p-p and Pb-Pb collisions at the LHC. The steeply falling 
shape above > p^ut^^cut _ g 2 jg reminiscent of that of parton fragmentation functions. 

In order to further show that similarity, the quark fragmentation functions into pions used in 
the calculation (either in the vacuum for p-p collisions or modified by the medium in Pb-Pb) 
at virtuality Q = p^"* are shown as solid lines. Note that at large 2;^^, no matching between 
correlation distributions and fragmentation functions is observed because of the onset of the 
fragmentation-photon channel. For more details, we refer the reader to [215]. 
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4 Outlook 

Parton fragmentation studies have been reviewed in this paper, trying to focus on the most 
recent advances of this active field. 

In the vacuum, no doubt that the (still preliminary) data from the B-factories and the more 
precise measurements to be carried out at RHIC with higher luminosities will prove essential to 
constrain further the fragmentation functions into light and heavy hadrons, specifically in the 
baryon sector. More refined experimental tests of the fragmentation processes at small- a; shall 
also be achieved at the Tevatron; equally impressive results can be expected by the ATLAS 
and CMS experiment at the LHC. On the theory side, a renewed interest has emerged recently 
with a systematic investigation of corrections beyond MLLA and for a variety of jet observables 
(single spectra, multiplicity distributions, 2-particlc correlations,. . . ). 

The understanding of QCD evolution in a medium is subject to a lot of attention. For years 
most of phenomenology has been based on the intuitive, yet rather simple, energy rescaling 
model. New interesting proposals have been suggested lately to go beyond that picture, with 
the modification of the evolution kernels in DGLAP or MLLA evolution equations. Despite those 
recent efforts, many questions are still open. It would be for instance interesting to investigate 
and understand how the QCD coherence effects observed for jets produced in the vacuum will 
manifest (or not) for jets traversing a dense and colourful medium; is there any angular ordering 
in QCD media? 

Finally, despite the now very precise RHIC data, it becomes increasingly clear that more 
differential observables will be needed to constrain the medium-modifications of fragmentation 
functions. Hopefully the possible full jet reconstruction in heavy-ion collisions and the measure- 
ments of photon-jet processes shall soon be performed and shed new light on our understanding 
of parton energy loss and the formation of QGP in heavy-ion collisions. 

It is a pleasure to thank the organizers of the Hard Probes 2008 conference, and in particular Nestor 
Armesto and Carlos Salgado, for having set up this exciting event. I am also indebted to David 
d'Enterria, Redamy Perez Ramos y Paloma Quiroga Arias for useful comments on the manuscript. 
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